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Abstract 

Previously, it has been shown that the direct correlation function for a Lennard-Jones fluid could 
be modeled by a sum of that for hard-spheres , a mean-field tail and a simple linear correction in the 
core region constructed so as to reproduce the (known) bulk equation of state of the fluid (Lutsko, 
JCP 127, 054701 (2007)). Here, this model is combined with ideas from Fundamental Measure 
Theory to construct a density functional theory for the free energy. The theory is shown to 
accurately describe a range of inhomogeneous conditions including the liquid-vapor interface, the 
fluid in contact with a hard wall and a fluid confined in a slit pore. The theory gives quantitatively 
accurate predictions for the surface tension, including its dependence on the potential cutoff. It 
also obeys two important exact conditions: that relating the direct correlation function to the 
functional derivative of the free energy with respect to density, and the wall theorem. 



1 



I. INTRODUCTION 



The insights of Van der Waals (VDW)^'^ still underlie much of the work on nonuniform 
fluids. One such insight is that the free energy of a fluid can be separated into two con- 
tributions: the first due to the short-ranged repulsion of all interaction models, the second 
due to a long-ranged attraction. The former is generally treated using an effect hard-sphere 
contribution and the latter by a mean-field approximation. This model is attractive be- 
cause of its simplicity and the prospect that it can be extended to complex systems such as 
molecular fluids, anisotropic interactions etc. However, while it can generally give reason- 
able qualitative predictions, such a simple model is seldom quantitatively accurate. It would 
therefore be highly desirable to modify the basic VDW model so as to give quantitatively 
accurate predictions if this can be done without compromising the basic simplicity of the 
model. The goal of this paper is to describe such a model and to show by application to 
the liquid-vapor interface and to fluids in confined geometries that it does indeed satisfy the 
twin requirements of accuracy and simplicity. 

In the language of density functional theory, the VDW model can be described as a 
particular approximation to the direct correlation function (DCF) of the bulk fluid: namely 
that it is the combination of a hard-sphere DCF and a mean-field tail. In a previous paper- 
it was noted that the failure of this model for the DCF lies primarily within the core region 
and that a good first order correction could be obtained by adding to the hard-sphere DCF a 
simple linear correction with coefficients adjusted to give continuity of the DCF and to give 
the correct thermodynamics. (The bulk thermodynamics are assumed to be known from e.g. 
thermodynamic perturbation theory or liquid state theory.) This is a satisfying approach as 
one of the motivations behind DFT (particularly the effective-liquid approaches) has been 
the idea that the DCF is a relatively simple function which could be easily approximated, in 
contrast to, say, the pair-distribution function which is highly structured. However, having a 
model DCF for the bulk fluid is only a first step towards describing non-uniform fluids: this 
model must somehow be used to construct a free-energy functional for nonuniform systems. 
If the DCF for an arbitrary nonuniform system were known, the free energy could be obtained 
immediately since the DCF is the second functional derivative of the free energy with respect 
to density. Thus one approach is to guess a generalization of the known DCF for a uniform 
system. In fact, since such functionals are known for the hard-sphere system and since the 
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mean-field tail is independent of density, the generalization need only be guessed for the 
core-correction which limits the problem. In ref.™, the usual ideas from DFT, involving the 
introduction of local densities into bulk expressions, were used to construct such a functional 
with mixed success: while the model gave reasonable predictions for the surface tension of 
the liquid-vapor interface, it suffered from some arbitrary elements due to the fact that 
the DCF and the proposed free energy functionals were internally inconsistent. That work 
therefore served to demonstrate the utility of trying to correct the DCF but did not fully 
address the problem of constructing a satisfactory free energy functional. 

In the case of hard-spheres, the problem of constructing a widely useful free energy 
functional has been solved in recent years by the development of Fundamental Measure 
Theory (FMT)^'^. Rosenfeld originally proposed FMT as a generalization of ideas from 
scaled particle theory«^, but for present purposes, one of the most interesting derivations 
of FMT is that of Kierlik and Rosinberg^. They get a theory essentially equivalent to 
Rosenfeld's by starting with the same general ansatz for the free energy functional (which, 
they point out, is the obvious generalization of the exact result known for one-dimensional 
systems) and requiring that the second functional derivative of this ansatz with respect to 
density give the Percus-Yevick DCF in the uniform limit. Thus, in this context, FMT can 
be seen as the result of a constructive exercise in which one starts with a known DCF, 
a particular ansatz for the free energy functional (which is exact in one dimension) and 
enforces the exact relation between the DCF and the free energy functional. Here, I propose 
that the the same constructive procedure be used to incorporate the core correction giving 
a fully consistent relation between the DCF and the free energy functional. Since the result 
is a straightforward modification of the hard-sphere contribution, while the mean-field tail 
is unchanged, the resulting model also satisfies the requirement for simplicity. This model 
will be referred to as the Modified Core VDW or MC-VDW free energy. 

The MC-VDW model requires as input the equation of state of the bulk fluid. This 
is in keeping with the view adopted in ref.^ that the main purpose of a model DFT is 
not to compute the properties of bulk fluids, which can be done very accurately using 
thermodynamic perturbation theory or liquid state theory, but to be used to calculate the 
properties of inhomogeneous fluids. This requirement is relatively mild compared to older 
DFTs that required knowledge of the direct correlation function of the bulk fluid as input. In 
the following, the Lennard- Jones system will be studied using both a perturabtive equation 
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of state and an empirical equation of state. 

In the next section, the details of the MC-VDW model are given. The third section 
consists of a comparison of the predictions of the model to data from computer simulations. 
The first comparison is the surface tension and density profile at the planar interface between 
coexisting liquid and vapor phases. It is shown that not only is the surface tension as a 
function of temperature accurately predicted, but so is the variation of the surface tension 
with the range of the potential. The second comparison is that of the structure of the fluid 
at a hard wall. It is noted that the theory satisfles the exact sum rule known as the Wall 
Theorem relating the density at the wall to the pressure far from the wall. The density as 
a function of distance to the wall is compared to the simulation results for several different 
temperatures. The flnal comparison is that of the structure of the fluid in slit pores. In all 
cases, the theory is found to be in good quantitative agreement with the simulations. The 
paper ends with a discussion of the results and of possible future developments. 



II. THEORY 

Given a collection of atoms interacting via a spherically symmetric pair potential, v{r) 
at flxed temperature T, chemical potential fi and external one-body fleld Vext (i"), the equi- 
librium density distribution p (r) is obtained by minimizing the functional 
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[p] = Jdv [p (r) Inp (r) - p (r) + /?/>< (r; [p]) - (3pp (r) + p (r) /^K.* (r)] (1) 



where j3 = l/ksT and is Boltzmann's constant^*^. The value of the functional at its 
minimum is the grand potential for the system. The only unknown term here is the excess 
free energy density, (r; [p]). It is related to the DCF for a nonuniform system via 

^^/^^^ 

c ri, rs; p = . ..... (2) 

6p (ri)5p(r2) 

where 

f3F^ = I f3r{T-[p])dr. (3) 

In ref.-, it was shown by a direct comparison to computer simulation that the DCF for the 
bulk fluid can be adequately approximated by a model of the form 

c (ri2; [p]) = CHS {ri2; d, [p]) + (^ao + Q {d - r^) - Pw (ri2) , (4) 



where the first term on the right is the hard-sphere DCF, d is the Barker- Henderson effective 
hard-sphere diameter, the constants and ai, which are functions of both density and 
temperature, are determined by requiring that this DCF gives the correct bulk free energy 
and that the DCF be continuous at the hard-sphere boundary ri2 = d (see ref.- and Appendix 
|X]for explicit expressions). The terms involving and Oi are referred to below as the "core 
correction". There are several reasonable choices for the tail function I3w{ri2) but here 
1 only consider the simplest choice, w {r-i2) = Q {r^ — d) v {r^). The idea is to use this 
model DCF for the homogeneous fluid together with the expression relating the DCF to the 
excess free energy functional to guide the construction of a functional that can be used for 
inhomogeneous fluids. 

In the simplest form of FMT-, the excess free energy for hard-spheres of diameter d is 
given by 

f^, (r; [p]) = s (r) $1 (r/ (r)) + $2 (v (r)) {s' (r) - (r)) + $3 (v (r)) s (r) {s' (r) - 3v' (r)) 

(5) 

where the explicit form of the algebraic functions are given in ref.-. The quantities 

ria (r) = (?7 (r) , s (r) , v (r)) are hnear functionals of the density 



j Wo(|r-ri|)p(ri)rfri. (6) 



The weights Wa (r) are 9(| — r), 5(| — r) and — r), respectively. In the uniform limit, 
the density becomes a constant p (r) ^ p and the quantity rjiy) — > TTpd^/6, which is the 
usual expression for the packing fraction. The other functionals become s (r) Trprf^and 
V (r) ^ 0. It is useful to introduce dimensionless quantities via 

*i (V) = {ri) (7) 

vr 



It is then straightforward to show that 
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{Gvh'l iv) + {6v? K iv) + (QV? h'i (r^)) - \x + ^x^) 
+ (2/i; (r/) + 4 (6r^) h', (r^) + 6 {Qr^f hi, (r/)) (i - ^x) 
+ (2/i2 (r7) + 6(6r/) hs{r])) (H 

(8) 



where x = ri2/d. For hard spheres, one way to determine the functions hi {rj) is to compare 
this expression to that for the Percus-Yevick DCF. 

Here, it is proposed to use the Modified-Core VDW free energy functional 

= j {(3f^s (r; [P]) + Pfcore (r; [p])) dr+^ J p (n) p (r,) 6 {r,, - d) v {r,,) dr.dr, (9) 
where the core correction is of the FMT form 

Pfcore (r; [p]) = ^Ji iv (r)) s (r)+^j2 iv (r)) {s' (r) - (r)) + ^j3 iv (r)) s (r) {s' (r) - 3v' (r)) 

(10) 

The functions ji (77) are introduced as analogs of the h-functions described above. They are 
determined by requiring that the model DCF be recovered in the uniform limit. Comparing 
eq.® and the core correction in eq.(jl]) gives 

= f; + (Qv) n + {Qvf Js (11) 

-ao -ai = 2j2 + 6 (677) jg 

-ao = j[ + 2iQ7])j', + 3i6vfji 

Its worth noting that use of these relations together with the explicit expressions for the 
coefficients ao and ai allows one to prove that these expressions do indeed reproduce the 
input bulk free energy in the uniform limit, 

l^PF^^re = l^PF^ - j^PF^HS -\-pje{r-d)v (r) dr 

= jiiv) + Qvj2{v) + (Qvfjaiv), (12) 

where N is the number of atoms, is the (known) bulk free energy of the fluid and F^g 
is the FMT hard-sphere free energy functional. The remainder of the derivation is given in 
Appendix |A] and only the flnal results will be given here. The function j^lrj) turns out to be 

J3 iv) = ^.[l {^l3F:ore - P^-pJ^f^Fcore ' j^PF^ore (0)) + 3^^" (v) " 3 J^' X^" iv) dv 

(13) 

where iv) = ^^^r^dCHS (^; P! d) is the hard-sphere DCF at the core boundary evaluated 
at the density corresponding to the packing fraction rj. (The hard-sphere DCF is completely 
determined by the hard-sphere FMT model: in the Rosenfeld model it is just the Percus- 
Yevick DCF^ .) Note that despite the fact that this expression results from the integration 
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of the differential equations in eq. (fTT|) . there are no integration constants. As shown in 
Appendix \^ any integration constants are forced to be zero by the requirement that the 
functions ji{ri) be finite at rj = 0. While this is not shown to be strictly necessary, it seems 
likely that divergences at zero density could lead to problems in inhomogeneous systems. 
Once jsirj) is determined, the remaining functions, ^1(77) and j2{v)j follow immediately using 
the second line of eq. (fTT!) and eq. (fT2|) . 

It is easy to see that, by construction, this free energy functional is consistent in the sense 
that eq.([2]), evaluated in the bulk limit, gives the assumed model DCF, eq.(jl]). 

III. COMPARISON TO SIMULATION 

In this Section, the results of the MC-VDW model will be compared to simulation results 
for the Lennard- Jones potential. 



which is typically used in Monte Carlo simulations. The only input required for the model 
is the equation of state of the bulk fluid. To test the idea behind the model, apart from 
other approximations, the very accurate, but empirical, 33-parameter equation of state of 
Johnson, ZoUweg and Gubbins (JZG)^ will be used. In order to illustrate the accuracy using 
more approximate methods that can be applied to other problems, the results using the first- 
order thermodynamic perturbation theory of Barker and Henderson (BH}^>^, and Weeks, 
Chandler and Anderson (WCA)^iiiii^ii^ will also be given. To provide some context. Fig. [1] 
shows the phase diagram calculated using all three of these together with simulation data. 
Clearly, the empirical equation of state is in close agreement with the simulations whereas 
the perturbative theories are reasonable at low temperatures but increasingly inaccurate as 
the critical point is approached, as is to be expected. 

Details concerning the numerical methods used in minimizing the free energy functional 
can be found in ref.~. In the following, temperature, T, and distance, z, will be expressed 
in reduced units as T* = T/e and z* = z/a respectively. A reduced density p* = pa^ 




(14) 



Some results will also be given for the truncated and shifted potential. 




(15) 
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will also be used. The hard-sphere contribution to the free energy was modeled using the 
White-Bear FMT functionaU^ii^ which is somewhat more accurate than the simplest FMT 
discussed above. However, for the class of problems considered here, it probably makes little 
difference which functional is used. 

A. The planar liquid-vapor interface 

The planar liquid-vapor interface is determined by minimizing the free energy functional 
for a value of the chemical potential corresponding to liquid-vapor coexistence and only 
allowing the density to vary in one direction (the z-direction) and with no external field. If 
the densities of the coexisting liquid and vapor are pi and p^,, respectively, then by definition 
of coexistence, the grand potentials of the two bulk phases are identical, [3VL (pi) = l3Q (p^). 
The excess free energy per unit surface, the surface tension, is then unambiguously defined 
as 

j = ^{Q[p]-Q (pj) (16) 

where A is the area of the surface perpendicular to the z-axis. The dimensionless surface 
tension is 7* = cr^7/e. Figure ([2]) shows the surface tension as a function of temperature as 
calculated from the theory and determined from simulations. Using the empirical equation of 
state, the calculated surface tension is consistent with the data given the scatter in the latter. 
The perturbative equations of state given reasonable values although the BH perturbative 
theory is somewhat superior to the WCA theory. Both are increasingly inaccurate at higher 
temperature due to their poor estimation of the critical point, where the surface tension 
vanishes. 

The figure also shows the calculated surface tension for a cutoff of r* = 2.5 which corre- 
sponds to that used in the simulations of Haye and Bruin^^. Using the empirical equation of 
state, the theory is somewhat less accurate in predicting the surface tension than in the case 
of the full potential, but the agreement is still reasonable. This is particularly the case when 
it is noted that the theory and simulation appear to extrapolate to slightly different critical 
temperatures which would account for most of the discrepancy and which indicates that it 
originates in the input equation of state. (Note that the modification of the JZG equation 
of state needed to account for the cutoff is not exact and is most inaccurate for very short 
cutoffs^.) The results using the approximate equations of state are similar to those found 
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with the full potential: the BH theory gives quantitatively better results, again probably 
due to the fact that it gives a better estimate of the critical point. 

The density profile at the interface is shown in Fig. ([3]) for several different temperatures 
and values of the cutoff. In all cases, the theory is in good agreement with the profiles 
determined from simulation!^. For T* = 0.7, near the triple point, and with a large cutoff, 
the theory predicts oscillations in the profile at the interface, as does for example the theory 
of Katsov and Weeks^^. However, the predicted oscillations here appear to be somewhat 
smaller than their prediction and more in line with the profiles observed in simulation^''. As 
the temperature is increased, the oscillations are quickly suppressed. A similar effect results 
from using a shorter cutoff. 

B. Hard Wall 

The next test is the determination of the density profile for a fluid in contact with a hard 
wall. In other studies, this comparison has been made using the simulation data of Balabanic 
et ai^. However, as this data is not readily available and as certain details such as the 
potential cutoff are unclear—, new simulations were performed using the Grand Canonical 
Monte Carlo method^^. The simulation procedure consisted of several steps. First, the 
chemical potential was estimated using the empirical equation of state. Then, an initially 
random configuration of atoms was run for 10^ attempted Monte Carlo moves with periodic 
boundary conditions. The simulation cell had length L in the x and y directions, and nL in 
the z direction where n is a parameter characterizing the geometry. This initial equilibration 
was followed by a further equilibration of 10^ attempted moves, but this time with hard walls 
at 2; = and z = nL. Finally, further runs of 10^ attempted moves were performed during 
which the density profile in the z direction was tabulated after every attempted moves 
using 200 equally-spaced bins where A^ is the expected average number of atoms. To control 
for the effect of system size, runs were performed at all densities involving approximately 
2000 atoms, corresponding to n = 2, and approximately 4000 atoms corresponding to n = 4. 
In all cases, the potential cutoff was r* = 4, corresponding (roughly) to that used in ref.— 
(see the discussion in ref.-°). 

Figures (jl]) and show a comparison of the theory (evaluated using the empirical 
equation of state) and simulation for temperature T* = 1.35 and chemical potentials corre- 
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spending to bulk densities of p* = 0.5, 0.65 and 0.82 respectively. (These conditions are the 
same as those used in previous studies^i22i2^). The theory clearly captures the very different 
qualitative behaviour across the range of densities and is quantitatively accurate for the 
lower densities. At the highest density, there is some difference in the calculated and mea- 
sured profiles away from the walls which can mostly be accounted for as an error in the phase 
of the oscillations with the theory. This is interesting as the FMT for hard-spheres gives, 
in the bulk, the Percus-Yevick structure and it has long been known that the Percus-Yevick 
approximation for the the pair-distribution funciton of bulk hard-spheres is also somewhat 
out of phase at high densities^. Given that the MC-VDW is a modified version of FMT, 
it is possible that the error in phase observed here is related. It should also be noted that 
this phase error would be partially or even totally negated if the points in the measured 
profile were incorrectly plotted at the left coordinate of the bins rather than at there center 
as has been done here. In all of the figures, it is clear that the agreement between simulation 
and theory is particularly good near the wall where the transition from drying at low bulk 
density to wetting at high density is correctly predicted. This is not an accident as it is 
known that DFTs of this form satisfy the exact sum rule p (0) = f3Pbuik where Pbuik is the 
pressure in the bulk far from the wall, a condition known as the wall theorem^ii^^. 

C. Slit Pore 

The final system considered is the Lennard- Jones fluid confined between two infinite walls 
(i.e. a slit pore). Unlike the previous case of a hard wall, the walls of the slit pore interact 
with the fluid via a modified Lennard- Jones potential intended to mimic the interaction 
between the fluid and a Lennard- Jones solid. The potential used here is the so-called 10-4-3 
potential of Steele, 



which is specifically meant to model the interaction between the fluid and a (100) plane 
in an FCC solid^^i^^i^^. Since the existing simulation data^'^'^ was obtained using very 
small numbers of atoms, new GCMC simulations with much larger systems were performed 
following the same protocol as for the hard wall. In the present case, the intermolecular 
potential was cutoff at r* = 6 and no cutoff was apphed to the wall potential. Figures ([MH]) 
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show a comparison between theory and simulation for the same conditions studied in ref.— 
and^, namely the chemical potential was set to the value corresponding to a bulk density 
of pI = 0.5925 and the temperature was T* = 1.2. The slit sizes in the five simulations are 
H* = 3, 4, 5, 6 and 10 and the simulations involved approximately 1300,2100,2400, 3000 and 
5000 atoms respectively. 

The calculations are in good agreement with the simulations. For the pores of intermedi- 
ate size, the calculations tend to overestimate the peak near the wall and to underestimate 
the amplitudes of the subsequent oscillations but are nevertheless reasonable. These re- 
sults are in broad agreement with previous models such as that of Tang and Wu^. Further 
calculations, not shown here, confirm the conclusions of Snook and van Megen^S that the 
profiles are insensitive to the value of the chemical potential so that errors in the equation of 
state are not as important as in the case of the liquid-vapor interface. Thus, the differences 
observed between theory and simulation must be attributable to the theory itself and can 
serve as a sensitive test for further improvements. 

IV. CONCLUSIONS 

In this paper, the modified-core VDW approximation for the direct correlation function 
has been used to construct a free energy functional based on ideas from Fundamental Mea- 
sure Theory. The resulting MC-VDW theory was shown to give accurate predictions for the 
surface tension of the liquid-vapor interface and the density profile near a wall and in slit 
pores. 

There are several advantages to the MC-VDW free energy functional. It has the practical 
advantage that it is no more complex than the mean field model constructed using the FMT 
for the hard-sphere contribution. However, unlike mean-field theory, it reproduces the input 
bulk free energy function so the bulk thermodynamics are automatically correct. It also 
satisfies two important exact relations. First, the exact relation between the free energy 
functional and the direct correlation function of the bulk phase is maintained: the second 
functional derivative of the free energy evaluated in the bulk phase gives the correct bulk 
DCF. Second, the wall theorem - relating the density at a hard wall to the pressure far from 
the wall - is satisfied. These relations are difficult to preserve in theories which are based 
on the introduction of the local density into bulk thermodynamic relations such as those 
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discussed in ref.-. The same comment applies to theories which attempt to locahze first 
order perturbation theory as discussed, e.g., in ref.—. In fact, the only theories based on 
these ideas which maintain the relation between the free energy functional and the bulk DCF 
are those which eliminate all density dependence beyond second order - which is basically 
the same as the earliest perturbative DFT of Ramakrishnan and Yussouff^^. 

As discussed in ref.-, the modified-core DCF was inspired by the work of Tang^^i^i^i^ 
on the First Order Mean-Field Approximation (FMSA). In that approach, the Ornstein- 
Zernicke equation is solved in a perturbative manner with the usual mean-field closure 
conditions. This gives an analytic result for the DCF of the bulk fluid which is similar to 
that used here, although the core correction is not linear in the spatial variable. It has the 
advantage that no input equation of state is needed, at least for Lennard- Jones interactions. 
However, the mean-field liquid state theory is known to be inaccurate in some circumstances^ 
and the fact that good results for the equation of state are obtained for Lennard- Jones by 
a first order solution seems somewhat fortuitous. Furthermore, it still remains to use the 
resulting bulk DCF to construct a DFT for inhomogeneous fluids. Different theories seem to 
be required, e.g., for confined fluids^ and for the liquid- vapor interface^. It is possible that 
the FMSA core correction could be used to construct a MC-VDW theory similar to that 
studied here, although the more complex analytic form of the core correction to the DCF 
will make this task more difficult. Nevertheless, this would seem a promising approach. 

The MC-VDW follows the philosophy described in the first paper in this series^: namely, 
that the bulk equation of state of fluids is well- understood using, e.g., thermodynamic per- 
turbation theory and that the real goal of DFT should be the construction of functionals 
for inhomogeneous systems. It would, of course, be more satisfying if, as in the case of the 
FMSA cited above, the equation of state could be derived from the theory. It is not impos- 
sible that this could be done using the present approach. Rather than using the equation 
of state to determine the functions defining the core correction (eqs jmfT3|) . these might be 
determined by imposing thermodynamic consistency between the resulting the free energy 
functional and, say, the internal energy calculated using the pair-distribution function given 
by the theory. This is the subject of ongoing work. 

A technical aspect of the MC-VDW model is that it is based on the simplest form of 
FMT as originally proposed by Rosenfeld^. This choice was made to avoid unneccessary 
complications in describing the model. However, in future work such as in application to 
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the sohd phase, it might prove necessary to use the more recent formulations of FMT— i^Si^, 
for the same reasons as occur in the case of hard-spheres^^. 

In summary, the present approach has the virtue of satisfying the relevant exact relations, 
of being conceptually and practically simple and of giving quite reasonable results for a 
variety of model problems. In particular, a single functional has been shown to be sufficient 
to describe both confined fiuids and the liquid- vapor interface. 
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APPENDIX A: DERIVATION OF MODEL EQUATIONS 

Using the second of eq. (fTT!) and eq. (fT3!) . one has that 

J2 iv) = -I (ao + ai) - 3 {Qr]) j, {r]) (Al) 

pV 

The third line of eq. (fTT|) becomes 

-do = j[ iv) + 2 M 32 iv) + 3 iv) (A2) 

= ^ (^/^^^r-e) - 3r7 (a', + a[) + 3 (ao + ai) - 72r]js (v) " SSr^'^Js (v) 



or, with some rearrangement 

o. , . 1 

di] 



(S6r]'j, {f]) - ^f3f:,re^ = 4ao + 3a, - 3r] {a', + a[) (A3) 



where /^^^g = ^/?-^core contribution of the core correction to the excess free energy 

density. 

For the chosen form of the tail-contribution to the DCF (which is independent of density), 
the explicit expressions for ao and a, are 

^ (4a„ + 3a0 = %<?> - + 4. H ,9» (r) r'dr (A4) 

3 dp^ dp' Jd 

ttQ + ai = —[3w {d) 
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where (p) = y(3F^ is the excess free energy density of the bulk fluid. Using these, the 
rhs can be written as 



4ao + 3ai — 3r] (oq + a[) (A5) 
Recognizing the contribution of the long-ranged tail to the bulk free energy density, 



fLi (^) Q{r-d) dr, 



(A6) 



this can be written as 



4ao + 3a, - 3r/ {a', + a[) = l^^-AIh±§ilL_J^ + 3v^chs {d-;p; d) (A7) 

Tia^ op at] 



or, since f^^^ = P - fns - ftau^ 



3 d-'Pl 



core 



d 

' di] 



(A8) 



d 

dr] 



Finally, integrating gives 



1^ 
2 dp 



1 

-PJcore 2 5p 



_f3fx _ J C01 

-PJcore 2 5p 



3r]—CHs {d-; p; d) 
dr] 



r d 

+ 3 r]—CHs{d-;p;d)dr] 
p=o Jo dr] 

(A9) 



where it is assumed that js (r^ = 0) = so as to avoid singulaties. 
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FIGURE CAPTIONS 



Fig. 1. (Color on line) The coexistence curve for the Lennard- Jones fluid as calculated 
using both the WCA perturbation theory, the BH theory and the empirical JZG equation 
of state. The full lines are the hquid-vapor coexistence curves, the dashed-lines are the 
spinodals and the symbols are the simulation data from ref.— (circles) and from ref.— . 

Fig. 2. (Color online)The surface tension as a function of temperature. The symbols 
are measurements from simulations (circles from ref.—, squares from TefM', diamonds from 
ref.— and triangles from^^). The lines are the results of the MC-VDW model evaluated with 
the JZG empirical equation of state (full line), the BH perturbation theory (dashed line) 
and the WCA theory (dash-dotted line). The lower curves and data are for a truncated and 
shifted potential with r* = 2.5 

Fig. 3. (Color online) Density profiles at the liquid- vapor interface calculated at different 
temperatures and values of the potential cutoff. From left to right, the curves correspond to 
T* = 0.7 and r* = 5.0, T* = 0.7 and r* = 2.5, T* = 0.8 and r* = 5.0, T* = 0.8 and r* = 2.5 
and T* = 1.1 and r* = 5.0. The symbols are the data reported in lefM and extracted from 
ref.— as the original is no longer available^. 

Fig. 4. (Color online) The structure of the fluid near a hard wall as determined from 
simulation (symbols) and the theory (lines). The simulations come from two runs each using 
cells with aspect ratio 1x1x2, circles, and 1x1x4, squares. The upper curve and data 
are for a chemical potential corresponding to bulk density p* = pa^ = 0.65 and the lower 
curve for density p* = 0.50. 

Fig. 5. (Color online) Same as fig. (jl]) except that the bulk density is p* = 0.85. 

Fig. 6. Comparison of the density distribution within slit pores of size H* = 3 and 
H* = 4 as calculated from the theory (lines) and as determined from simulation (symbols). 

Fig. 7. Same as fig. for H* = 5 and H* = 6. 

Fig. 8. Same as fig. (El) for H* = 10. 
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